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Abstract 

Nonlinear dynamics of a bouncing ball moving vertically in a gravi- 
tational field and colliding with a moving limiter is considered and the 
Poincare map, describing evolution from an impact to the next impact, 
is described. Displacement of the table is approximated in one period by 
four cubic polynomials. Results obtained for this model are used to elu- 
cidate dynamics of the standard model of bouncing ball with sinusoidal 
motion of the limiter. 



1 Introduction 

In the present paper we study dynamics of a small ball moving vertically in a 
gravitational field and impacting with a periodically moving limiter (a table). 
This model belongs to the field of nonsmooth and nonlinear dynamical systems 
[JJ [21 121 H] ■ In such systems nonstandard bifurcations such as border-collisions 
and grazing impacts leading often to complex chaotic motions are typically 
present. It is important that nonsmooth systems have many applications in 
technology El CD HI IS] • 

Impacting systems studied in the literature can be divided into three main 
classes: bouncing ball models [101 El [H], impacting oscillators [T3] and im- 
pacting pendulums [TH H] , see also [I] . In dynamics with impacts it is usually 
difficult or even impossible to solve nonlinear equation for an instant of the next 
impact. For example, in the bouncing ball models the table's motion has been 
usually assumed in sinusoidal form, cf. [12] and references therein. This choice 
of the limiter's motion leads indeed to nontractable nonlinear equation for time 
of the next impact. To tackle this problem we proposed a sequence of models in 
which periodic motion of the table is assumed (in one period of limiter's motion) 
as a low-order polynomial of time |15j . It is thus possible to approximate the 
sinusoidal motion of the table more and more exactly and conduct analytical 
computations. Carrying out this plan we have studied several such models with 
linear, quadratic and cubic polynomials [TBI [T71 [TBI IT9] . 



1 



In the present work we conduct analytical and numerical investigations of 
the model in which sinusoidal displacement of the table is approximated in one 
period by four cubic polynomials. We shall refer to this model as Mc- Si- 
multaneously, we study the standard dynamics of bouncing ball with sinusoidal 
motion of the limiter, referred to as Ms- We hope that rigorous results ob- 
tained for the model Mc cast light on dynamics of Ms- It should be stressed 
that results obtained for the model Ms can be compared with experimental 
studies, see [3UJ US H2 for the early papers, summarized in [53], and [23] for 
recent work. 

The paper is organized as follows. In Section 2 a one dimensional dynamics 
of a ball moving in a gravitational field and colliding with a table is reviewed 
and the corresponding Poincare map is constructed and models of the limitcr's 
motion Mc and Ms are defined. Bifurcation diagrams are computed for Mc 
and Ms- In Sections 3, 4 and 5 a combination of analytical and numerical 
approach is used to investigate selected problems of dynamics in models Mc 
and Ms- More exactly, fixed points and their stability are discussed in Section 
3, birth of low velocity n-cycles is investigated in Section 4 and birth of high 
velocity 3-cycles is studied in Section 5. In Section 6 the case of N impacts in 
one interval of the limiter's motion is studied for the model Mc- We summarize 
our results in the last Section. 

2 Bouncing ball: a simple motion of the table 

Let a ball moves vertically in a constant gravitational field and collides with a 
periodically moving table. We treat the ball as a material point and assume 
that the limiter's mass is so large that its motion is not affected at impacts. 
Dynamics of the ball from an impact to the next impact can be described by 
the following Poincare map in nondimensional form [25] (see also Ref. [TT] 
where analogous map was derived earlier and Ref. |12) for generalizations of the 
bouncing ball model): 



where Ti denotes time of the i-th impact and Vi is the corresponding post-impact 
velocity while Aj+i = Tj+i — Tj. The parameters 7, R arc a nondimensional 
acceleration and the coefficient of restitution, < R < 1 [5], respectively and 
the function Y (T) represents the limiter's motion. The limiter's motion has 
been typically assumed in sinusoidal form, Ys(T) = sin(27rT). Equations (TJJ 
and Y — Ys lead to the model Ms- This choice of limiter's motion leads to 
serious difficulties in solving the first of Eqns.([T]) for Ti + i, thus making analyt- 
ical investigations of dynamics hardly possible. Accordingly, we have decided 
to simplify the limiter's periodic motion to make dial) solvable. The function 



V i+ i 



7 Y(T i )-A? +1 + A; +1 ^, 

-RVi + 2RA i+1 + 7 (1 + R) Y (Ti+i) , 



(la) 
(lb) 
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Yc (T): 



Yc(T) 




< T < { 

1 < f < - 

4-^2 

1 < T < - 

2 — 1 ^ 4 

| < T< 1 



(2) 



/ 2 (T) 

/a(r) 



= (32tt - 128) T 3 + (-16tt + 48) T 2 H 
= (128- 32tt)T 3 + (-144 + 32tt)T 2 
= (128 - 32tt) f 3 + (-240 + 64tt) f 2 



2ttT (3a) 
f (48 - IOtt) f - 4 + 7T (3b) 
f (144 - 42tt) f - 28 + 9tt (3c) 



(32tt - 128) T 3 + (336 - 80?r) T 2 + (-288 + 66tt) T + 80 - 18tt (3d) 



approximates I5 = sin(27rT) on the intervals [k, k + 1], fe = 0, 1, . . . , with T = 
T — [T\ , where [a^J is the floor function - the largest integer less than or equal 
to x. The model Mc consists of equations (JTJ) , @, © with control parameters 

R, 7. We shall also need velocities of the limiter, defined as gi (T) = (T), 
i = 1 4. 




Figure 1: Bifurcation diagram for the model M.c, R = 0.85. 

In Fig. 1 above we show the bifurcation diagram with impact times (blue) 
and velocities (black) versus 7 computed for growing 7 and R — 0.85. It follows 
that dynamical system Mc has several attractors: two fixed points which after 
one period doubling give rise to chaotic bands and two other fixed points which 
go to chaos via period doubling scenario. There are also several small attractors. 



3 



We shall investigate some of these attractors in the next Section combining 
analytical and numerical approach (general analytical conditions for birth of 
new modes of motion were given in [26 ). 

We show below the corresponding bifurcation diagram for the sinusoidal 
motion. 
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Figure 2: Bifurcation diagram for the model M.s, R = 0.85. 



Similarity of Figs. 1, 2 suggests that analytical results obtained for model 
A4c shed light on the problem of sinusoidal motion, Ms- 



3 Fixed points and their stability 

We shall first study periodic solutions of the model Aic with one impact per k 
periods and T £ (0, ^) since it is suggested by the bifurcation diagram that 
they are stable. Such states have to fulfill the following conditions: 



V n+1 = V n = Vj k/1) , T n+1 =T n + k = ri fc/1) + k (k = 1, 2, . . .) , (4) 



where: 



X) e (0, i) , K (fc/1) > lY Cl (Ti k ^) . (5) 
The demanded (stable) solution is given by 

(fe/i) 



5 = ef^-^iyVV-e) 2 +65(^-4), (i? = ii T |)(6a) 
K (fc/1) - k. (6b) 



*(s) 
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Since T» £ [0, 1] we demand that T* > and it follows from (|6ap that physical 
solution appears for lower critical value 7 > 7c?/ 1 ' where: 

(fe/l) _ _L!^R (7) 
lcn,C — 2ir l+R- \') 

We have checked by stability analysis that the solution (I5a|) . (|6bl) is stable 

for 7 > 7^j c"; i- e - when it is physically acceptable. To determine upper critical 
value of 7 when dynamics looses stability we put into ((TJ) : 



Ti = T§M + e u T i+1 = T§™ + k + e l+1 , (8) 
Vi = V, + ^ = k + Hi, V l+ i = V, + i-h+i = k + fi l+1 , (9) 

with Y (T) given by ([3]), and keep only terms linear in perturbations £j, £i+\, 
fii, /Ltj+i of the fixed point to get: 



(Z\)-{l(l + R)9i(T*) 1 ^$&™ -r)(*) (10) 

where T* = T§£\ /1 (T) is given by © and 5l (T) = ^/j (T). 
Since the characteristic polynomial is: 

X 2 + aX + [3 = 



a - 4^4 (tt - 6) 2 + 6fc (tt - 4) ^gJT (1 + i?) 2 7 - i?2 _ 1 (n) 
/3 = i? 2 

application of the Shur-Cohn criterion ([27 ): 

(3 < 1 (12) 
\a\ < [3 + 1 

leads finally to the localization of the fixed points ([6]), 7^ q < 7 < with: 



(A/1) _ 6fe( 7 r-4)(fl 2 -l) + V36fc2( 7 r-4) 2 (l-^ ) 2 +4(7r _ 6) 2 (1+jR2) 2 
Tcr 2 ,C - 8(^=6^(1+1? ' ^ J 

In Fig. 3 stability regions in (i2, 7) plane for the .Mc model arc shown. In 
the case of the model Ais we have: 



'cr 2 ,S — 'cr 2 ,C> K 1 ^) 
(fe/1) _ V^^(l-^) 2 +4(l + fl 2 ) 2 

Tcr 2 ,s - 2¥^T+77P ' • K<1 > Uoj 

see [IS] (note that in [25] we used Y(T) = sin(T) rather than Y S (T) = sin(2vrT) 
and it follows that all values of the control parameter A must be rescaled, 7 = 
, A , 2 ) and stability regions are very similar to those of model A4c, cf. Fig. 4. 

{27T) 
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Figure 3: Stability regions in the (R, 7) plane, model Mc- 
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Figure 4: Stability regions in the (-R, 7) plane, model Ms- 
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4 Birth of low velocity k - cycles 



In this Subsection we shall study birth of low velocity k - cycles which can be 
seen in the bifurcation diagrams, Figs. 1, 2, for 7 > 0.03 and V < 1. In the 
case of such cycles Ti,T 2 , • • • , Tfc e (0, 1) and Tk+\ — 1 = T\, Of course, it is 
possible to follow periodic orbits backwards, i.e. iterating the map ([1]) until the 
convergence to the k - cycle is achieved for some initial condition and some 7. 
Then the value of 7 is decreased (slightly) and the map is iterated again (until 
convergence is obtained) with the previously computed fc - cycle as the initial 
condition. This method, although leads to determination of the critical value 
of 7 at which the k - cycle disappears for decreasing 7 (or is born for growing 
7) but is time-consuming and not very effective due to very poor convergence 
near the threshold. 

On the other hand, analytical conditions for birth of k - cycles are found 
below. In what follows theorems about differentiation of implicit functions [28j 
will turn out useful since Eqn. (|laj) defines Tj_|_i implicitly. Consider equation: 



F(T 1 ,T 2 ) = 0, (16) 

which defines dependence of, say, T 2 on T±, see |28j where necessary and suf- 
ficient assumptions are given. Then it follows from implicit function theorem 
that: 

dT 2 F{ 



where F[ = J£, F> 2 = 

In a more complicated case, equations: 

F(T 1 ,T 2> r 3 ) = 0, G (Ti,T 2 ,T 3 ) = 0, (18) 

define T 2 and T3 as functions of T\ under appropriate assumptions. We can now 
compute derivatives with respect to T\ as [28] : 



•Ml F i 5 I det I ^ G 2 



dT 2 _ _ \ F 3 G' 3 J aTs _ V F l G"l 



9Tl " 7 Ft GL \ ' aTl " / Fi G 



(19) 



with F[ = ^jr, F 2 = F3 = jjjr and analogous notation for i = 1, 2, 3. 

4.1 Low velocity 2 - cycle in the model Ale 

Numerical tests show that a 2 - cycle fulfilling conditions T\ € (0, j) , T 2 <E 
(|, f ) and T 3 = T x + 1 is stable. This 2 - cycle can be seen in the bifurcation 
diagram in Fig. 1 for 7 > 0.0366 and V x = 0.51, V 2 = 0.55 (R = 0.85). 
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Equations to determine T\, T 2 and Vi, V 2 are shown below: 

7/3(12) = 7/1 (Ti) - (T 2 - Ti) 2 + (Tfe - Ti) Vi 

V 2 = -Wi + 2i?(T 2 -Ti)+7(l + i?) 5 3(T 2 ) 

7/i(T3-l) = 7/3(T 2 )-(T 3 -T 2 ) 2 + (T 3 -T 2 )y 2 

V 3 = -RV 2 + 2R{T 3 -T 2 )+ 1 (l + R)g 1 (T 3 -l) 

T 3 = T x + l 

V 3 = Vi 



(20) 



where /, (T)'s and <7i (T)'s are defined in Eqn. ([3]) and the text below. 

We were able to simplify Eqns. (|20[) signihcantly obtaining equation for 
A = T 2 - Ti only: 

F ( A ) = E- = o^ AJ = ' (21) 

where c?i 's are given in the Appendix. Numerical computations suggest that the 

(2) (2) 
2 - cycle appears for 7 = % r C and fixed i?, where % r C is a critical value, as 

_____ 

a double (and stable) solution of Eqns. (j2"0")l . For 7 > % r C there are two real 
solutions, one stable (seen in the bifurcation diagram) and another unstable. 

(2) 

On the other hand, for 7 < 7V. c the solutions are complex conjugated and thus 

i 2 \ 

unphysical. Moreover, at 7 = % r ( j the stability matrix has unit eigenvalue. 
Therefore this is a tangent (saddle-node) bifurcation, see 29 for elementary 
discussion of the tangent bifurcation in the logistic map when the 3 - cycle is 
born. All other cycles discussed in our paper are also born in tangent bifurcation. 

To determine critical value of the parameter 7 let us note that double solution 
of the polynomial equation (|2"Tj) is also the solution of G (A) = where G (A) = 
■j^F (A). For example, solving for R = 0.85 the system of equations: 

F(A) = ]T 9 =0 ^ AJ=0 ' (22) 
G (A) = E- = i^ iAi " 1 = °' (23) 

we get 7^ c = 0.036 617 052 682 892 250 62, A cr = 0.634 279 960 677 747 355 95 
(and many other, unphysical solutions) in perfect agreement with numerical 
computations, see also Fig. 1. 

Alternatively, we can use implicit function theorem. Solving the second and 
fourth equations in (|20|) for V\ , V 2 we get 



„ 7 (1 + R) (R93 (T 2 ) - g! (Ti)) + 2R (1 + R) (r 2 - Tj) - 2R 

Vl = ^TT^ 

7(1 + R) Cggi (Tj) - g 3 (T 2 )) - 2R (1 + R) (T 2 - 7j) + 2i? 2 
2 -1 + i? 2 

and 

F (Ti, T 2 ) £ 1 f 1 (Ti) - 7/3 (T 2 ) - (T 2 - Tx) 2 + (T 2 - Ti) Vi = 

G (Ti, T 2 ) I 7 / 3 (T 2 ) - 7/j (Ti) - (Ti + 1 - T 2 ) 2 + (Ti + 1 - T 2 ) V 2 = 



(24) 



(25) 
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We can, in principle, solve the equation F (Ti, T 2 ) = to determine T 2 (Ti) 
and demand that ^yr-G (Ti,T 2 (Ti)) = to obtain condition for double root: 

F{T U T 2 ) - 

G{T U T 2 ) = (26) 
^-G(Ti,T 2 ) = g§r-G(Ti,T 2 ) + -^jrG{Ti,T 2 ) fpf- = 

where the derivative ^£ is computed from Eqn. (fT7|) . Eqns. (|26|) provide 
analytical condition for the onset of the 2 - cycle. They are too complicated 
to be solved analytically but can be solved numerically for a fixed value of R 
or 7. For example, for R = 0.85 we compute the critical value of 7 and the 
critical 2 - cycle: T x = 8. 167 748 882 344 294 132 7 x 1CT 2 , T 2 = 0.715 957449 
501 190 297 28, j^} c = 3. 661 705 268 289 225 062 x 10~ 2 in perfect agreement 
with solution of Eqns. (|22"]l. ([25]). 

4.2 Low velocity 2 - cycle, model .Ms 

We can apply this result to the case of sinusoidal motion. First of all, there is 
analogous 2 - cycle with T-y € (0, |) , T 2 6 (|, |), which appears at — 
0.0346, see Fig. 2. It can be thus assumed that this 2 - cycle is also born as a 
double solution. The corresponding equations of the 2 - cycle are of form: 



7sin(2^T 2 ) = 7sin(27tTi)-(ra-Ti) 3 + (r a -ri)7i 

V 2 = -RV 1 +2R(T 2 -T 1 )+~f(l + R)2Trcos(2irT : 



2) 



7 sin (2tiTi) = 7 sin (2ttT 2 ) - (T t + 1 - T 2 f + (Tj + 1 - T 2 ) V 2 



2 



(27) 



Vi = -RV 2 + 2R(T 1 + 1-T 2 )+j(1 + R)2ttcos{2ttT 1 ) 
Solving second and fourth equations of (|2"T|) for V\, V 2 we get: 



Vi = [RTi ~ R T 2 + 77T cos (27rTi) - jRtt cos (2ttT 2 ) + ^ 



(28) 



V2 = iffl ( --RT! + i?T 2 - 7i?7rcos (2tiTi) + 77rcos (27rT 2 ) - ^ 
The problem is thus reduced to the system of two equations for T\, T 2 

F(T U T 2 ) = 7sin(27rTi)-7sin(27rT 2 )- A 2 + yiA = 

G{Tx,T 2 )= 7sin(27rT 2 )-7sin(27rTi)- A 2 + y 2 A = ( 29 ) 
A = T 2 -T l! A = Ti + l-T 2 



We couldn't solve the system of equations (|28|) . (|29|) analytically. Analytical 
condition for double root of these equations, i.e. for the beginning of the 2 - 
cycle, are again provided by Eqns. (|26|) . (fl7|) with functions F, G defined in 
(|29|) . Solving now these equations numerically for R = 0.85 we get the critical 
value of the parameter 7 and values of dynamical variables of the critical 2 - 
cycle: 7^ = 3. 458 072 636 337 462 017 6 x 10~ 2 , T x = 7.171 236 860 717 641 004 
3 x 10~ 2 ' T 2 = 0.706 981 761 358 463 856 03. 
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(2) 

Numerical computations show that at 7 = 7L, s there is indeed a double 



lcr,S 



solution of (|2~T1) . two real solutions for 7 > % r s (one stable, another unstable) 

and complex solutions for 7 < %r . These considerations describe and explain 
the birth of the corresponding 2 - cycles. 



4.3 Low velocity 3 - cycle in the model Aic 

We have found numerically that a 3 - cycle satisfying conditions 7\ £ (0, j) , T 2 e 
(|, §J , T 3 g (|, l) and T4 = T\ + 1 is stable. This attractor is seen in the 
bifurcation diagram near the 2 - cycle for 7 > 0.0452, R — 0.85, cf. Fig. 1. The 
3 - cycle variables fulfill equations: 

7/3(72) = 7/i(7 1 i)-(T 2 -T 1 ) 2 + (T 2 -T 1 )y 1 

V 2 = -RV l +2R(T 2 -T 1 )+j(l + R)g 3 (T 2 ) 
7/473) = lf 3 (T 2 )-(T 3 -T 2 ) 2 + (T 3 -T 2 )V 2 

V 3 = -RV 2 + 2R{T 3 -T 2 )+ 1 {l + R)g i {T 3 ) 
7 /i(T 4 -l) = 7/4 (T 3 ) - (T 4 - T 3 ) 2 + (T 4 - T 3 ) V3 

V4 = -RV 3 + 2R (T 4 - T 3 ) + 7 (1 + i2) 51 (74 - 1) 

T 4 = + 1 

^4 = Vi 

Equations (j3"0|) can be simplified. We can solve the second, fourth and sixth 
equations for Vi, V 2 , V 3 to get 



(30) 



Vi = rgi (Ti) + TR 2 g 3 (T 2 ) - TR 9i (T 3 ) + aT x + bT 2 - cT 3 + d 
V 2 = -TRgi (Ti) + Tg 3 (T 2 ) + Ti? 2 54 (T 3 ) - cT x + aT 2 + bT 3 - dR 
V 3 = TR 2 9l (7i) - TRg 3 (T 2 ) + Tg 4 (T 3 ) + bT x - cT 2 + aT 3 + dR 2 

r _ 7 _ 2fl(l-fl) , _ 2fl 2 2fl . _ 2fl 

1 ~~ fi 2 -i?+l> " ~~ fl 2 -_R+l' " — fl 2 -_R+l' L — ft 2 -,R+l' u— 1+7? 3 



(31) 



The problem is thus reduced to three equations for impact times Ti, T 2 , T 3 only: 



F (Ti, T 2> T 3 ) i 7 /i (Ti) - 7/3 (T 2 ) - A 2 + A1V1 = 

G (T x , T 2 , T 3 ) i 7/3 (T 2 ) - 7 / 4 (T 3 ) - A 2 + A 2 y 2 - 

if (T U T 2 , T 3 ) i 7/4 (T 3 ) - 7/1 (71) - A 2 + A3V3 = 

A 1 =T 2 -T 1 , A 2 =T 3 - T 2 , A 3 =T 1 + 1 

where V\, V 2 , V 3 are known functions of impact times, cf. (j3Tj) . We were unable 
to solve Eqns. (|30|) analytically. However, it is possible to write down condition 
for the onset of the 3 - cycle since it follows from numerical computations that 
the 3 - cycle is born as a double root of Eqns. (|30|). The condition for the 
double root is -£rH (21, T 2 (2i) ,T 3 (T x )) = and hence the condition for the 



(32) 

-T 3 
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onset of the 3 - cycle is: 



F (Ti, T 2 ,T 3 ) = 

G(T U T 2 ,T 3 ) = (33) 



H(T U T 2 ,T 3 ) = 



where the derivatives , are computed from JT2J). Solving these equa- 
tions numerically for R = 0.85 we get critical value of the control parameter 
7 and the critical 3 - cycle: ^ c = 4. 518 834447 846 807 553 9 x 10" 2 , T x = 
0.103 931 597 153 962 754 97, T 2 = 0.635 200 266 830 198 212 15, T 3 = 0.848 760 321 
631572 414 99. 



4.4 Low velocity 3 - cycle in the model M.$ 

We can apply this result to the case of sinusoidal motion. First of all, there is 
analogous 3 - cycle with T x £ (0, |) , T 2 £ (±, §) , T 3 £ (§, l) and T 4 = T x + l, 

which appears at Y cr — 0.04499, see Fig. 2. We can thus expect that this 3 
- cycle is also born as a double solution. Dynamical variables of the 3 - cycle 
obey equations: 



7 sin(2^T 2 ) = 7 sm(27rTi) - (T 2 - T^ 2 + {T 2 -T 1 )V 1 

V 2 = -RV 1 +2R(T 2 -T 1 )+ 1 (l + R)2ircos(2irT 2 ) 

7 sin(2^r 3 ) = 7sm(2^T 2 )-(T 3 -T 2 ) 2 + (T 3 -T 2 )V 2 

V 3 = -RV 2 + 2R(T 3 -T 2 )+j(1 + R)2ttcos(2ttT 3 ) 

7 sin (27rTx) = 7 sin (2^T 3 ) - (T x + 1 - T 3 f + (T x + 1 - T 3 ) V 3 

Vi = -RV 3 + 2R{T 1 + \-T 3 )+ 1 {1 + R)2ttcos{2-kT 1 ) 

Solving second, fourth and sixth equations for V\, V 2 , V 3 we get 

Vi = a (R 2 C 2 - RC 3 + d)+b (R (T 2 - T x ) +T 1 -T 3 ) + c 
V 2 =a (R 2 C 3 - Rd +C 2 )+b (R (T 3 ~ T 2 ) +T 2 - T{) - cR 
V 3 = a (R 2 C 1 - RC 2 +C 3 )+b (R (T x - T 3 ) +T 3 - T 2 ) + cR 2 

a = 7 1 + R 3 ' ' ^ = — 1+R 3 ' C = 1 + R 3 ' ^ = C0S ^ 7r -^i (* = 1j 2, 3) 

and we have to solve equations for impact times only: 

F (Ti,T 2 , T 3 ) i 7 sin (2tiTi) - 7 sin (27rT 2 ) - A 2 + AiVi = 

G (Ti, T 2) T 3 ) i 7 sin (2ttT 2 ) - 7 sin (2^T 3 ) - A 2 + A 2 F 2 = 

H (Ti, T 2 , T 3 ) I 7 sin (2ttT 3 ) - 7 sin (2ttTi) - A 2 + A 3 V 3 = 
Ai = T 2 - T u A 2 = T 3 - T 2 , A 3 ^T 1 + 1-T 3 



(34) 



(35) 



(36) 



Equations (|35|) . (|36|) are too complicated to be solved analytically. However, 
it is possible to write down condition for the beginning of the 3 - cycle since it 
follows from numerical computations that the 3 - cycle is born as a double root 
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of Eqns. (fM|) . More exactly, we have to solve Eqns. (f3"3")l for functions F, G, H 
denned in ©. We thus get for R = 0.85 the critical value 7^^ = 4. 498 669 496 
445 746 754 8 x 1CT 2 and the critical 3 - cycle, Ti = 9.514 258132 574445 543 
3 x 1CT 2 , T 2 = 0.633 092 075 481 873 314 56, T 3 = 0.848 082 849 264 211 982 09. 



5 Birth of high velocity 3 - cycles 

High velocity 3 - cycles are very characteristic of bouncing ball dynamics. They 
accompany all fixed points and are seen in the bifurcation diagrams around 
V = 1, 2, 3, ... , see Figs. 1, 2. In the case of such cycles li,T 2 ,T 3 G (0, 1) 
and T4 — k = T\ . 



1+R» ' 1+R 3 ' l+Tf 3 



and 



5.1 Model M c , V ^ 1 

We start with such 3 - cycle with V = 1 which appears in the model M.q for 
7 > 0.042, see Fig. 1 with impact times 7\ G (±, §) , T 2 - 1 G (0, |) , T 3 - 1 G 
(0, j) , T4 — 1 = T\. The corresponding equations are: 

7/1 (T 2 - 1) = 7/2 (Ti) ~ (T 2 - T1) 2 + (T 2 - Ti) Vi 

V2 = -fl7i + 2R (T a - Ti) + 7 (1 + i2) 51 (T 2 - 1) 

7/i(T 3 -l)= 7/i(r 2 -l)-(T 3 -T 2 + l) 2 + (T 3 -T 2 + l)V 2 , , 

V 3 = -W 2 + 2#(T 3 -T 2 + l)+7(l + fl)i ?1 (T 3 -l) ^ ' 

7/2 (Ti) = 7/1 (T 3 - 1) - (Ti + 2 - T 3 ) 2 + (Tr + 2 - T 3 ) V, 

Vi = -RV a + 2R (Tt + 1 - T 3 + 1) + 7 (1 + i?) <? 2 (Ti) 

Solving equations for Vi, V 2 , V 3 we get 

Vi = a lg2 (Ti) - ai?7.gi (T 3 - 1) + ai? 2 7 .gi (T 2 - 1) - 2aRA 1 + b 

V 2 = ja 9l (T 2 - 1) - ai? 752 (Ti) + a 7 i? 2 5l (T 3 - 1) - 2aJL4 2 - 6 

l/ 3 -a7.9i(T 3 -l)-ai?7 5l (T 2 -l) + a7i? 2 52 (T 1 ) + 2ai?A 3 + c 

Ai = -RT 2 + RTi - Ti + T 3 , A 2 = -RT 3 + RT 2 -T 2 +T 1 ( 38 ) 

A 3 = -RT a + RTi -T 2 + T 3 

_ (l+R) 1. _ 2R(2-R) _ 2R(l+2R 2 ^ 

(X -, . . " '' 



^i (Ti,T 2 , T 3 ) * 7 / 2 (Ti) - 7 A (T 2 - 1) - A 2 + A x V l = 

^2 (Ti, T 2 , T 3 ) * 7/x (T 2 - 1) - 7 /i (T 3 - 1) - A 2 + A 2 V 2 = (3g) 

^3 (Ti, T 2 , T 3 ) I T /x (T 3 - 1) - 7/2 (Ti) - A 2 + A 3 V 3 = 
Ai = T 2 - Ti, A 2 = T 3 - T 2 + 1, A 3 = Ti + 2 - T 3 

Analytical condition for the onset of this 3 - cycle is given by Eqns. (f3"3"| with 
F, G, H given by (|39|) . Solving these equations numerically for R = 0.85 we get 



jg* 1 ' = 4. 184 258 672 013 445 046 3x 10~ 2 and Ti = 0.292 944 344 346 867 579 94, 
T 2 = 1.113 429 439 708 681876 6, T 3 = 1.179198 676 090 233 556 9. For R cr . c = 



12 



0.685 101 194 and 7 cr ,c = 0.056 81 9 493 there is smooth transition to the state 
Ti G (\, \) , T 2 € (§', l) , T 3 - 1 G (0, ±) , T 4 - 1 = T v For i? > i? cr , c this 
transition occurs for 7 > 7 cr; c- 



5.2 Model M s , V ^ 1 

In the case of sinusoidal motion described by the model Aic we can see the 3 - 
cycle with V = 1 in Fig. 2 for 7 > 0.045 with impact times Ti G (0, |), T 2 -1 G 
(i' I)' -^3 — ^ e (0) z)) ^4 — 1 = ^X- Dynamical equations read: 

7 sin (2ttT 2 ) = 7 sin (2tiTi) - (T 2 - Ti) 2 + (T 2 - Ti) Vi 

F 2 = -RVi + 2R (T a - Ti) + 7 (1 + ii) 2tt cos (27rT 2 ) 



7 sin (27tT 3 ) = 7 sin (27rT 2 ) - (T 3 - T 2 + l) 2 + (T 3 - T 2 + 1) F 2 

V3 = - W 2 + 2i? (T a - T 2 + 1) + 7 (1 + iZ) 2tt cos (2ttT 3 ) 

7 sin (2ttTi) = 7 sin (27rT 3 ) - (Ti + 2 - T 3 f + (Ti + 2 - T 3 ) F 3 

Pi = - RV 3 + 2i? (Ti + 2 - T 3 ) + 7 (1 + iZ) 2tt cos (2ttTi) 



and 



2 77 r(l+fl) , _ 2R(1+R) _ 2R(R-2) , _ 2R.(l+2R 2 ) 
1 + R :i ' — 1 + ' C — l+R 3 ' fl — 1 + i? 3 



F (Ti, T 2) T 3 ) i 7 sin (2tiTi) - 7 sin (2ttT 2 ) - A 2 + AiF = 

G (Ti, T a , T a ) i 7 sin (2ttT 2 ) - 7 sin (2^T 3 ) - A 2 + A 2 V 2 = 

# (Ti, T 2 , T 3 ) i 7 sin (2ttT 3 ) - 7 sin (2ttTi) - A§ + A 3 V 3 = 
Ar = T 2 - Ti, A 2 = T 3 - T 2 + 1, A 3 = Ti - T 3 + 2 



(40) 



Solving equations for V\, V2, V 3 we get 

Vi = aR 2 cos 2vrT 2 - aR cos 27rT 3 + a cos 2ttTi - k4i - c 

y 2 = ai? 2 cos 2?rT 3 - aR cos 2ttTi + a cos 27rT 2 - bA 2 + ci? 

y 3 = aR 2 cos 2ttTi - aR cos 27rT 2 + a cos 27rT 3 + bA 3 + d 

M = -RT 2 + RTi - Ti + T 3 , A 2 = -iZT 3 + i?T 2 - T 2 + Ti ( 41 ) 

A 3 = JZTi - i?T 3 + T 3 -T 2 



(42) 



Conditions for the onset of 3 - cycle are given by Eqns. (13"3l with functions 
F, G, H defined in (|42]) . Solving these equations for i? = 0.85 we obtain: 
7^' 1} = 4.514 020 805 615 479 8341 x 10~ 2 and Ti = 7.439 906 099 929 247 941 
lx 10~ 2 , T 2 = 1.154 226 438 813 261 052 6, T 3 = 1.357478 350 324075 728 9. For 
Rcr,s — 0.691 964 922 5 and 7 CJ -,s = 0.055 974 756 there is smooth transition to 
the state Ti G (±, |) , T 2 G (§, l) , T 3 - 1 G (0, ±) , T 4 - 1 = T\. For R > R cr , s 
this transition occurs for 7 > 7 cr ,s- 



6 N impacts in one period of limiter's motion 
and chattering in the model A4c 

In the bouncing ball dynamics chattering and chaotic dynamics arise typically, 
see |30[ 131] where chattering mechanism was studied numerically for sinusoidal 
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motion of the table. Due to simplicity of our model analytical computations are 
possible. 

We shall consider a possible course of events after grazing. 



6.1 First interval: T„ T l+l e (0, J) 

Let Ti, Tj_|_i G (0, j). In this case we get from Eqns. (fTaj) . (|3ei|) A 4 +i = and: 

A (±) _ h°i ( T ») + 1 (, + L 128 7 (4- 7tTW \ 
1+1 647(4-.) ^ 1/ (i 7aim) + l)V' (43) 
at (T) = **f x (T) , 

and A| + J is the solution describing chattering (obviously, Wi must be small 
enough so that expression under the square root be non- negative) . The de- 
nominator in (j43l) can be written as ^701 (Ti) + 1 = —967 (4 — n) (Ti — T$ 
where: 

(1) -l + 16 7 (-3 + 7r) 

96 7 (-4 + 7r) ' 1 ' 

and we check that T$ < 2 occurs for 7 > 7 cr' ) = 0.043 731. Therefore for 
7 < 7cr the grazing ball will stay on the table forever. Let us assume that 
a ball sticks to the table for some time T g < < \. At critical point 

Ti = Tcr^ and Vi = 751 (Ti), equations ((Taj) , (|3a|) have the degenerate triple 
solution Ti + i — Tcr . For 7 > jir and T > T^ the solution aL_ \ is no longer 
valid since ^70.1 (Tj) + 1 < and A\ + \ < what is physically unacceptable. The 
solution A^ is also unacceptable and thus the ball has to jump to another time 
interval, (j, 5), (5, |) or further. We shall now consider the first possibility. 

Let us assume that the ball grazes at Tj = T^ and thus its velocity is that of 
the table, Vi = 751 (Ti). We thus have to solve equation for the jump: 



7/2 (T i+ x) = 7/1 (Ti) - (T 4+ i - T 4 ) 2 + (T i+1 - Tj) V t 
Ti = G (0, i) , T l+1 G (i, |), V, = 7.91 (T) . 

The solution of (gSJ) is 



(45) 



T (^2) _ (^^)(4 + 2, + ^j + l 1 + 2 f + ^2 

-4 + 7T r 96(-4 + 7r) 7 ' V ; 

It follows that the interval (t£\ ri 1 ^ 2 ^ is the forbidden zone. The solution 

(|50"j) is valid for 7 < 7^^ 2) = 0.057 102 since for 7 > ji? we have T^ 1 ^^ > 
~ contradicting assumptions. For 7 > 7 cr >2 ^ we thus have to consider the 
following equation for the jump to time interval (h, |): 

7/3 (T i+1 ) = 7 /i (T,-) - (T i+ i - T^) 2 + (T !+1 - Tj) V5, 
Ti - T C W G (0, i) , T l+1 G (|, |), V, = 75i (31) . 
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Solution of Eqn. (l47l) . X — Tcr >3 \ fulfills the following cubic equation: 

a X 3 + (b + 617) X 2 + (c + C17 + c 2 7 2 ) X + d + di7 + ^7 2 + 47 3 = 
a = 32, 6 = 1, h = -240 + 64tt, c = ^ , ci = 1 - ±7T 
c 2 = 304tt- 552 - d = 57^48, dx = -576 7T+ 

d 2 = - 1^ + i d 3 = fr^ 3 - ^tt 2 + 364tt - 444 

x " (48) 

with Tcr >3 ' ) = 7 (_4+ tt) • Eqn. (|52|) has acceptable solutions, i.e. such that 

T c ( r^ 3) € (£, |), for 7 < 7 ^^ 3) = 0.087 308 825. 
6.2 Second interval: 2}, T m € (J, 5) 

We have to consider now chattering in the interval |). Let us thus consider 
that Ti,T i+ i 6 Q, |). It follows from equations (fTa|) . ((3b| that A i+ i = and 



A (±) _ |7«2 (Ti) + 1 / t ± ( T ') + 1 A , 128 7 (4- 7r)Wi 



1+1 64 7 (4-tt) ^ |i 7 a 2 (T,) + l|y (§ 7 a 2 (T*) + l)7 ' 
a 2 (T) = (T) . 

(49) 

The solution describing chattering is A^tJ for ^702 (Ti) + 1 > 0. The de- 
nominator can be written in form ^702 (Ti) + l = 967 (4 — tt) (T — X^ 2 -* ) where: 



_ l l-»4 7 + 32 7 » 
96 7 (7T — 4) 

It follows that chattering is thus possible for Ti € {Tcr\ j>J ■ 
6.3 Third interval: T;,T i+ i e (|, |) 

Let us suppose now that Tj,Tj+i e (i, |). It follows from equations flTa|) . (j3cj) 
that A,+i = and 



a (±) _ 17^3 (T«) + 1 / t ± f 7°3 ( r ») + 1 /, , 128 7 (4- 7r)Wi 



4+1 64 7 (4-.) ^ |i 7a3 (T,) + l|y faMTJ + l)*)' 
as (T) = jg*f a (T) ■ 

(51) 

The solution describing chattering is A^tJ for ^703 (Ti) + 1 > 0. The de- 
nominator can be written in form ^703 (Ti) + 1 = 967 (4 — n) ^T — Tc?^j with: 

1 -2407 + 64 7 7r + 1 , 

T -=m 7( ,-4) <0 - 5 ( ^ >0) < (52) 

and it follows that Tc? cannot belong to (i, |) interval for positive 7. Chat- 
tering is thus possible in the whole interval (i, |) since ^703 (Ti) + 1 > 0. 
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6.4 Fourth interval: T it T i+1 G (|, l) 

Let us assume finally that two subsequent impacts occur in the last quarter- 
period and Ti, Tj+i 6 (|, l). In this case the solution A i+ i = of equation 
(fTaj) is always present and this equation can be easily solved. We thus get from 
Eqns. (Taj, (H3| A. i+ i = and: 



(±) ^4(^ + 1 / / 128 7 (4-7r)^ 



04 (T,) - ^ACO. (54) 

where 704 (Tj) is the acceleration of the table, 04 (T) = -^sA (T) and Wi is a 
relative velocity, Wi = V\ — 7^-/4 (Tj). In the case of chattering the appropriate 
solution is A[ + \ since for Wi — ► we have A^ , \ — > 0. The denominator can 
be written as ±704 (T 4 ) + 1 = -96 7 (4 - tt) (t - Tj?) where 



m 1 167(21 - 5?r) + 1 
r £ ) = 96 7(4-^) ^ (7>0) - (55) 



We do not have to worry that the denominator in (|53p may vanish since 
the condition ^7a 4 (Tj) + 1 = cannot be fulfilled for Tj G (|, l) and 7 > 0. 
Therefore chattering is possible in the whole interval (|, |) since ^704 (T) + l > 
0. 

6.5 Grazing: a homoclinic orbit 

Let us assume that the ball grazes at Tj = with velocity Vj = jgi (Ti) (i.e. 
it has velocity of the table) and that the value of 7 is such that the ball jumps. 

Let us next assume that the ball grazes at some time T* in the interval 0, T$ . 

For growing 7 it will happen eventually at 7 = 7, that T* = • Then for 
larger values of 7 the ball after grazing and jumping returns with chattering into 
the ^0, Tcr^ interval but it will not graze, i.e. V* will be larger than 7171 (T»). 

We have computed numerically the critical value as 7* = 0.058 348 6. Therefore 
for 7 > 7* long transients can be expected after grazing. 

In Fig. [5] below bifurcation diagram is shown with initial conditions on the 
grazing manifold. For 7 < 7* the grazing manifold is globally attractive. Indeed, 
in this parameter range the bifurcation diagram is empty (we show attractors 
different than the grazing manifold only). On the other hand, for 7 > 7, the ball 
jumps at T = Tj* and then either grazes eventually or settles on some attractor 
after a long transient. Just after the threshold there is a very irregular, probably 
chaotic attractor. 
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Figure 5: Bifurcation diagram for the model Aic, R = 0.85. Initial conditions 
are on the grazing manifold. 



7 Summary 

We have studied dynamics of a bouncing ball impacting with a periodically 
moving limiter within two frameworks of the table motion: Aic and Ms defined 
in Section 2. Stability conditions of fixed points have been determined and 
results for the models Aic and Ais have been compared. Then we have found 
that low- velocity fc-cycles as well as high-velocity 3-cycles are generically born 
in tangent bifurcations. Moreover, we have been able using implicit functions 
theorems, to write down conditions for the onset of these cycles and solve them 
numerically. Analytical conditions for the onset of such cycles are new. 

Finally, the case of N impacts in one interval of the limiter 's motion has 
been studied within the Aic model. Equations for N impacts in one period 
of limiter's motion were found and simplified significantly, making analysis of 
chattering and grazing possible. We have found, combining analytical and nu- 
merical approach, the grazing homoclinic orbit which appears at 7 = 7* and 
gives rise to a very irregular, probably chaotic attractor, cf. Fig. [5] We expect 
that analogous attractor exists in the model Ais- 
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A Equation for the 2 - cycle 

In the Appendix coefficients of the polynomial (|2Tj) are listed. 

d 9 = 2 16 3 3 7 3 (-4 + tt) 4 (1 + Rf 

d s = -2 12 3 3 7 2 (1 + Rf (-4 + tt) 3 ((72tt - 288) 7 + 1) 
d 7 = 2 6 3 3 7 (1 + R) 4 (-4 + tt) 2 (a 2 R 2 + ai R + a ) 

a 2 = (-660487T+ 82567r 2 + 132096) 7 2 + (256tt- 1024)7+1 
a x = (-138 2407r+ 17 2807r 2 + 276 480)7 2 + 2 + (5127r- 2048)7 
a a = (-660487T+ 82567r 2 + 132096) 7 2 + (256tt- 1024)7+1 

d 6 = -8(1 + Rf (-4 + tt) (b 2 R 2 + hR + b ) 

b 2 = (-19 160 064tt 2 + 76 640 256tt+ 1596 672tt 3 - 102 187 008) 7 + 
(-3024 + 756tt) 7 + (1361 664 - 679 680tt + 84 864tt 2 ) 7 2 + 1 

61 = (-241 532 928 + 3773 9527r 3 + 181 149 6967r- 45 287 4247r 2 )7 3 + 
(183 552?r 2 - 1469 952tt + 2944 512) 7 2 + (1512tt - 6048) 7 + 2 

b = (-19 160 0647T 2 + 76 640 2567r+ 1596 6727r 3 - 102187 008)7 3 + 
(1361 664 - 679 680tt + 84 864tt 2 ) 7 2 + (-3024 + 756tt) 7 + 1 

d b = 8 (1 + Rf {c 2 R 2 + Cl R + c ) 

c 2 = (l80 983 808- 180 983 8087r- 11 319 5527r 3 + 707 9687r 4 + 67 885 0567r 2 ) 7 3 + 

(2951 424tt - 3953 664 + 61 056tt 3 - 734 976tt 2 ) 7 2 + 

(17 928 + 11187T 2 - 8952tt) 7 + 3tt - 12 
ci = (2202 3687r 4 + 565 014 528 - 564 682 7527r+ 211636 2247r 2 - 35 254 2727r 3 )7 3 + 

(163 584tt 3 - 10 561 536 - 1967 616tt 2 + 7893 504tt) 7 2 + 

(- 16 176tt + 32 400 + 2020tt 2 ) 7 - 24 + 6tt 
c = (706 4327r 4 + 68 106 2407r 2 - 181979 1367r- 11327 2327r 3 + 182 310 912)7 3 + 

(2951 424tt - 3953 664 + 61 056tt 3 - 734 976tt 2 ) 7 2 + 

(17 928 + 1118tt 2 - 8952?r) 7 + 3tt - 12 
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d 4 = -4(1 + R) 2 (e 4 R 4 + e 3 R 3 + e 2 R 2 + e x R + e ) 

e 4 = (50 264 064 - 31818247r 3 + 201 3447r 4 - 50 264 0647r+ 18 929 6647r 2 )7 3 + 

(1854 720tt - 2515 968 - 456 960tt 2 + 37 632tt 3 ) 7 2 + 

(29 160 - 14 520tt + 1810tt 2 ) 7 - 36 + 9tt 
e 3 = (-29 354 496tt 3 + 473 112 576 + 176 357 376tt 2 + 1834 496tt 4 - 471 453 696tt) 7 3 + 

(-2864 640tt 2 + 11 566 080tt - 15 593 472 + 236 928tt 3 ) 7 2 + 

(6160tt 2 + 99 360 - 49 440tt) 7 + 30tt - 120 
e 2 = (-68 8967r 4 + 20 570 1127r+ 1156 8007r 3 - 21689 856- 7310 8807r 2 ) 7 3 + 

(-27 039 744 - 4981 248tt 2 + 412 416tt 3 + 20 086 272tt) 7 2 + 

(140 400 - 69 840tt + 8700tt 2 ) 7 + 42tt - 168 
ei = (486 383 616+ 1819 1 367r 4 + 178 569 2167r 2 - 29 431 2967r 3 - 481406 9767r)7 3 + 

(-2864 640tt 2 + 11 566 080tt - 15 593 472 + 236 928tt 3 ) 7 2 + 

(6160tt 2 + 99 360 - 49 440tt) 7 + 30tt - 120 
e = (-3220 2247r 3 + 193 6647r 4 - 55 240 7047r+ 20 035 5847r 2 + 56 899 584) 7 3 + 

(1854 720tt - 2515 968 - 456 960tt 2 + 37 632tt 3 ) 7 2 + 

(29 160 - 14 520tt + 1810tt 2 ) 7 - 36 + 9tt 

d 3 = -16 (1 + Rf (f 4 R 4 + f 3 R 3 + f 2 R 2 + hR + f ) 

U = (3905 280- 233 088tt 3 + 14 OOOtt 4 - 3891 456tt+ 1437 2647r 2 ) 7 3 + 

(-16 128tt + 27 648 + 3072tt 2 - 192tt 3 ) 7 2 + 

(-3132 + 1548tt - 192tt 2 ) 7 + 8 - 2tt 
h = (340 032tt 3 - 6027 264- 2066 112tt 2 - 21 504tt 4 + 5709 312tt)7 3 + 

(115 968tt 2 - 479 232tt + 663 552 - 9408tt 3 ) 7 2 + 

(-622tt 2 - 10 152 + 5016tt) 7 - 5tt + 20 
f 2 = (-68 8967r 4 + 20 570 1127r+ 1156 8007r 3 - 21689 856 - 7310 8807r 2 )7 3 + 

(1714 176 + 308 736tt 2 - 25 344tt 3 - 1257 984tt) 7 2 + 

(-12 312 + 6072tt - 752tt 2 ) 7 - 6tt + 24 
/1 = (-9234432- 17 792tt 4 - 2600 640tt 2 + 358 592tt 3 + 8114 6887t) 7 3 + 

(115 968tt 2 - 479 232tt + 663 552 - 9408tt 3 ) 7 2 + 

(-622tt 2 - 10 152 + 5016tt) 7 - 5tt + 20 
fa = + (-225 088tt 3 + 15 600tt 4 - 2854 656tt+ 1206 864tt 2 + 2522 880) 7 3 + 

(-16 128tt + 27 648 + 3072tt 2 - 192tt 3 ) 7 2 + 

(-3132 + 1548tt - 192tt 2 ) 7 + 8 - 2tt 
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d 2 = 2 (1 + Rf (g 4 R 4 + g 3 R 3 + g 2 R 2 + 9l R + g ) 

34 = (9372 672 - 534 528tt 3 + 31 552tt 4 - 9206 784tt + 3349 440tt 2 ) 7 3 + 

(-152 064 + 2656tt 3 - 31 200tt 2 + 120 384tt) 7 2 + 

(-3456 + 1632tt - 196tt 2 ) 7 + 36 - 9tt 
,g 3 = (-1544 4487r 3 + 24 55 1 424 + 9416 4487r 2 + 93 1847r 4 - 25 049 0887r)7 3 + 

(-114 432tt 2 + 419 328tt - 497664 + 10 176tt 3 ) 7 2 + 

(-11927T 2 - 19 872 + 9696tt) 7 - 12tt + 48 
32 = (I33 2487r 4 - 25 214 9767r- 1969 9207r 3 + 21 731328+ 10 696 3207r 2 )7 3 + 

(1838 592 + 318 336tt 2 - 25 728tt 3 - 1321 344tt) 7 2 + 

(-12 096 + 5760tt - 696tt 2 ) 7 - 18tt + 72 
.gi = (l2 607 488+ 107 0087r 4 + 7425 7927r 2 - 1475 3287r 3 - 16 091 1367r)7 3 + 

(-114 432tt 2 + 419 328tt - 497 664 + 10 176tt 3 ) 7 2 + 

(-11927T 2 - 19 872 + 9696tt) 7 - 12tt + 48 
30 = (-515 328tt 3 + 35 392tt 4 - 6718 464tt+ 2796 480tt 2 + 6054 912) 7 3 + 

(113 472tt - 138 240 - 30 624tt 2 + 2720tt 3 ) 7 2 + 

(-3456 + 1632tt - 196tt 2 ) 7 36 - 9tt 

di = 2 (1 + Rf (h 4 R 4 + h 3 R 3 + h 2 R 2 + hxR + ho) 

hi = (-570247r 3 + 38727r 4 - 746 4967r + 3113287r 2 + 663 552) 7 3 + 

(-13 824 + 416tt 3 - 4128tt 2 + 13 248tt) 7 2 + 

(360tt - 48tt 2 - 648) 7 + 3tt - 12 
h 3 = (148 224tt 3 - 3151 872- 1021 824tt 2 - 7552tt 4 + 2985 984tt) 7 3 + 

(-216 576tt - 4800tt 3 + 276 480 + 56 064tt 2 ) 7 2 + 

(152tt 2 - 1248tt + 2592) 7 
h 2 = (-25 920tt 4 + 5474 304tt + 395 136tt 3 - 4976 640- 2223 9367T 2 ) 7 3 + 

(-32 640tt 2 + 132 480tt + 2688tt 3 - 179 712) 7 2 + 

(720tt - 96tt 2 - 1296) 7 + 6tt - 24 
h x = (-497 664- 10 6247r 4 - 579 4567r 2 + 132 8647r 3 + 995 3287r) 7 3 + 

(-216 576tt - 4800tt 3 + 276 480 + 56 064tt 2 ) 7 2 

(152tt 2 - 1248tt + 2592) 7 
ho = (-57 0247r 3 + 38727r 4 - 746 4967r + 3113287r 2 + 663 552) 7 3 + 

(20 160tt + 352tt 3 - 27648 - 4704tt 2 ) 7 2 + 

(360tt - 48tt 2 - 648) 7 + 3tt - 12 



20 



d = k 6 R 6 + k 5 R 5 + k 4 R 4 + k 3 R 3 + k 2 R 2 + k x R + k 

k 6 = (-235 8727r 2 + 608 2567r- 25287r 4 - 580 608 + 401287r 3 ) 7 3 + 

(4272tt 2 + 17 280 - 400tt 3 - 14 976tt) 7 2 + 

(216 - 120tt + 16tt 2 ) 7 - tt + 4 
fc 5 = (-7744tt 4 + 1 14 048tt 3 - 1327 104+ 1492 992tt- 622 656tt 2 ) 7 3 + 

(-32tt 3 - 6912 + 3456tt - 288tt 2 ) 7 2 + 

(-240tt + 32tt 2 + 432) 7 
fc 4 = (-74567r 4 + 995 3287r- 746 496 - 478 3687r 2 + 98 8807r 3 ) 7 3 + 

(18 432tt - 24 192 + 368tt 3 - 4560tt 2 ) 7 2 + 

(-360tt + 648 + 48tt 2 ) 7 + 12 - 3tt 
k 3 = (552 960tt - 256 896tt 2 + 52 480tt 3 - 442 368 - 3968tt 4 ) 7 3 + 

(-480tt + 864 + 64tt 2 ) 7 
fc 2 = (-588 960tt 2 + 1492 992tt+ 102 7207r 3 - 66887r 4 - 1410 048) 7 3 + 

(-42727T 2 + 14 976tt + 400tt 3 - 17 280) 7 2 + 

(-360tt + 648 + 48tt 2 ) 7 - 3tt + 12 
fci = (-7744tt 4 + 1 14 048tt 3 - 1327 104+ 1492 992tt- 622 656tt 2 )7 3 + 

(288tt 2 - 3456tt + 32tt 3 + 6912) 7 2 + 

(-240tt + 32tt 2 + 432) 7 
fc = (38 848tt 3 + 442 368tt - 199 008tt 2 - 359 424 - 2784tt 4 ) 7 3 + 

(-368tt 3 + 4560tt 2 - 18 432tt + 24 192) 7 2 + 

(216 - 120tt + 16tt 2 ) 7 + 4 - tt 
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